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Y* ■ Abstract 

The Euler-Poincare approach to complex fluids is used to derive multiscale equations for 

£^ computationally modelling Euler flows as a basis for modelling turbulence. The model is 

O ,■ based on a kinematic sweeping ansatz (KSA) which assumes that the mean fluid flow serves 

as a Lagrangian frame of motion for the fluctuation dynamics. Thus, we regard the motion of 

a fluid parcel on the computationally resolvable length scales as a moving Lagrange coordinate 

for the fluctuating (zero- mean) motion of fluid parcels at the unresolved scales. Even in the 

simplest 2-scale version on which we concentrate here, the contributions of the fluctuating 

l/~) \ motion under the KSA to the mean motion yields a system of equations that extends known 

results and appears to be suitable for modelling nonlinear backscatter (energy transfer from 

CO . smaller to larger scales) in turbulence using multiscale methods. 
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1 Multiscale approaches 

1.1 Dealing with microstructure dynamics in turbulence 

The history and present state of the art of multiscale approaches to fluid turbulence modelling are 
recounted in a number of excellent sources, including [EE03, Hou2005, PS2008, EH2009, Tr2010]. 
For the present purpose of modelling multiscale fluid turbulence, we mention Multiscale Finite 
Elements Methods [EH2009] and the Heterogeneous Multiscale Method [HMM]. Both are general 
methodologies for numerical computation of problems with multiple scales. The methods rely on 
an efficient coupling between the different macroscopic and microscopic physical models. The key 
to the efficiency of such an approach is the possibility that the microscale model may not need to 
be solved over the entire computational domain, but only over small selected regions near where 
data estimation is carried out. Examples of applications include complex fluids, micro-fluidics, 
solids, interface problems, stochastic problems, and statistically self-similar problems. 

The present work applies the standard multiscale method to Euler's fluid equations, then 
combines the results with ideas from the geometric mechanics of complex fluids, in order to create 
a new two-scale model of ideal incompressible flow. This is accomplished by: (i) making a slow- 
fast spatial decomposition of the fluid velocity; (ii) performing Lagrangian averaging in Hamilton's 
principle using this decomposition; and then (iii) applying the Taylor hypothesis in assuming that 
fluctuations are convected by the mean flow. The last step treats the flow trajectories of the mean 
flow as Lagrangian coordinates for the dynamics of the fluid fluctuations. The result is a two-scale 
extension of the Lagrangian-averaged Euler alpha equations of [HMR1998] which was the basis for 
the Lagrangian-averaged Navier-Stokes alpha (LANS-a) equations of [FHT2001, FHT2002]. 



1.2 Convection of microstructure 

Homogenization techniques were applied in [MPP1985] to obtain an averaged equation for the 
large scale features of highly-oscillatory solutions of the three-dimensional incompressible Euler or 
Navier-Stokes equations. The following initial value problem was treated: 

d t u + (u ■ V)u = — Vp, with V ■ u = 

and with highly-oscillatory initial data (e<l) 

u{x,0) =U(x) + W (x,-) . 

Multiscale expansions were constructed for both the velocity field and the pressure, under the im- 
portant assumption that the microstructure is adverted by the mean flow. Under this assumption, 
the following multiscale expansion for the velocity field was constructed: 

_, . , . (6{x,t) t x \ /6(x,t) t x \ ^ / o\ 

u e {x, t) = u(x, t)+wl -^- L - L , -, -, x, t\ + em ( - ±J - L i -, -, x, t J + O (e 2 ) . 

The pressure field p e was expanded similarly. This form of the solutions for the velocity and 
pressure fields were shown to be consistent with the fluctuation quantity 6 being advected by the 
mean velocity as a Lagrangian coordinate. Namely, 

d t 9 + u-V9 = 0, 9(x,0) = x. (1) 

The additional vector variable 9 is the back-to-labels map, or inverse map for the three-dimensional 
incompressible Euler equations. Modelling the effects of the rapid small scales on the slower 
large-scale solutions of the three-dimensional Euler and Navier-Stokes equations constitutes the 
closure problem in turbulence theory. The work in [MPP1985] provided some understanding of 
the interactions of the small scales with the large scales and it characterized the back-to-labels 
map as a form of fluctuating microstructure attached to the mean flow. The closure problem for 
turbulence was not solved, however, because the solutions for the functions u and w of the rapidly 
oscillating variables turned out not to be unique. The uniqueness problem was addressed by 
imposing additional assumed constraints that led to large-scale averaged equations that resembled 
the then-popular k-e closure model of turbulence. 

Convection of microstructure of the two- and three-dimensional incompressible Euler equations 
has also been studied from a related but different viewpoint from [MPP1985] in a series of recent 
papers, culminating in [HYR2008]. These papers used a new approach that does not assume 
that the fluctuations follow the characteristics of the mean flow. This new approach enabled the 
derivation of a well-posed cell problem which was used to obtain an effective large-scale average 
equation. However, the numerical solution still required solving for the back-to-labels map at each 
time step. 

The present work builds on the results of this previous research, and obtains systems of mul- 
tiscale equations that are purely Eulerian and avoid the necessity of solving for the back-to-labels 
map. The approach is based on the Euler-Poincare theory of complex fluids developed during the 
past decade, see e.g. [HT2009] and references therein, also in a series of other papers [GBTRPH]. 
The complex fluids approach is based on geometric mechanics and on the application of reduction 
by symmetry to Hamilton's principle for ideal continua [HMR1998]. The corresponding Hamil- 
tonian approach to complex fluids was developed earlier by Dzyaloshinskii and Volovik in their 
famous paper [DzVol980]. Later, a Hamilton's principle approach was followed to establish the 
variational approach to various soft matter systems in [GBTRPH], where various complex fluids 
were studied in different contexts, including dissipative dynamics. 



Plan of the paper and main results 

fi) Section 2 discusses the two-scale fluid flow decomposition that we use and introduces Taylor's 
hypothesis as a kinematic sweeping ansatz (KSA). 
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Section 3 formulates the two-scale model by applying geometric methods for Hamilton's 
principle that were borrowed from complex fluids theory. 

Section 4 discusses conservation laws for total subgrid scale circulation and helicity. 

Section 5 develops resolved-scale models based on advection laws for subgrid scale quantities. 
These models are similar to the stretched subgrid scale vortex model of [MP1997]. 

Section 6 discusses a subgrid scale fluid model with feedback between the two scales. 

Section 7 contains a summary and conclusions of the paper. 

There are also five Appendices that explain details of proofs and derivations sketched in the 
body of the text. 



2 Flow decomposition and Kinematic Sweeping Ansatz 

2.1 Decomposition of periodic vector functions 

Lemma 1 (Decomposition Lemma [HYR2008]) 

Any periodic function u(x) : R 3 — > 1R 3 that admits a Fourier series expansion on the unit cube 
[0, l] 3 may be decomposed into a sum of two periodic functions 

u(x) = ui(xi) + u 2 (xi,x 2 ) with xi,x 2 G[0, l] 3 (2) 

in which u 2 has zero mean in x 2 G [0, l] 3 . That is, J u 2 (xi, x 2 ) d 3 x 2 = 0. 

The proof follows from manipulation of Fourier series and is provided in the present notation 
in Appendix A.l. We will refer to ui and u 2 in equation (2) as the velocities of the mean flow 
and the fluctuations, respectively. The xi dependence in u 2 (xi,x 2 ) may be regarded as the slowly 
varying envelope of the rapidly fluctuating component of the velocity. The fluid pressure has a 
similar decomposition. 



Potential for applications in numerics. We will use Lemma 1 to express Euler's fluid equa- 
tions in terms of two scales (xi,x 2 ) = (x, x/e), in which we will regard e G [0,1] as a fixed 
parameter. In a computational setting, for example, the value of e could be assigned by the choice 
of grid size used to resolve the large-scale solution Wi(x), while regarding the remainder w 2 (x, x/e) 
as the unresolved, subgrid-scale part of the solution. Thus, e = 1/N in this case corresponds to 
the fractional length scale that one would like to resolve on a computational grid. The limit e — > 
would pick up more and more small-scale components of the solution, and it 2 (x, x/e) would tend to 
zero in the limit. That is, letting e — > resolves more and more small scales into the solution. In 
that limit, the computation would become a direct numerical simulation which attempts to resolve 
all scales. However, that limit is not our objective. 



The challenge for us here is to introduce a computable fluid model that describes the effect of 
subgrid scales on the resolvable scales xi at a fixed nonzero value of e. Our approach will be to use 
the hypothesis of convected fluid micro structure [MPP1985, HYR2008]. For this, we will regard 
the motion xi(£) of a fluid parcel on the resolvable length scale as a moving Lagrange coordinate 
for the motion x 2 (xi(i), t) of fluid parcels at the unresolved scales. We will apply this assumption 
by using methods of geometric mechanics. After some kinematic considerations for convection of 
fluid microstructure, the result will be an Eulerian multiscale fluid model. 

2.2 Lagrangian averaging: fluctuations are swept by the mean 

The main hypothesis of the present paper is that fluctuations are swept by the mean flow, but 
they are not slaved to the mean flow as in Large Eddy Simulation (LES) turbulence modelling. 
In particular, we assume that the fluctuations are swept by the large-scale motion and they have 
zero mean in the Lagrangian frame moving with the large-scale velocity. 

Our interpretation of fluctuations is that they undergo their own evolution, which in turn 
depends on the Lagrangian fluid parcel traveling with the mean flow. That is, at each mean 
Lagrangian label if)(x.\ ), there corresponds another Lagrangian label f7 Xl (x;> ) that is associated 
to fluctuation dynamics. More specifically, we introduce the following definition. 



KSA: Fluctuations are swept by the mean. That is, while the Lagrangian label x^ for the mean 



Definition 1 (KSA: Kinematic Sweeping Ansatz) 

KSA: Fluctuations are swept by the mean. 1 
flow is taken to its current position by the map 

xi = v(x! 0) ) , 

the fluctuation label x 2 transforms according to 

The subscript notation in the above formula emphasizes that the Lagrangian map r) for the fluctu- 
ations depends on the Eulerian position x x of the mean flow parcel. The maps i/> and r] also both 
depend explicitly on time. Moreover, both the mean flow and fluctuations are assumed to undergo 
incompressible dynamics, so the maps if) and r\ each preserve their respective volume elements 

d 3 xS 0) = d 3 x x and d 3 x 2 0) = d 3 x 2 . 

As a result of the KSA, the velocity of a fluctuation as it is swept along a Lagrangian mean 
trajectory Xi(t) may be written in the form 

^ = |„ xi (x<°») + u 1 .V 1 ^,(xf), (3) 

where u^X!) = ip(tp~ (xi)) is the Eulerian mean velocity and Vi stands for d/dx^. 

Remark 1 The main difference between the present two-scale sweeping approach and other ap- 
proaches in the literature lies in the second material term ui(xi) ■ Vi77 Xi , which shows that the 
fluctuations are dragged by the mean flow in a Lagrangian sense, so that Lagrangian mean trajec- 
tories Xx(t) become Lagrangian variables for the fluctuation dynamics. 



The fluctuations will be constrained to have zero mean. However, the question arises of how 
the mean should be taken. It is clear that the precise quantity possessing zero mean is the 
Eulerian fluctuation velocity u 2 (x 1 ,x 2 ,£). However, one must ask in which sense the latter has 
constant zero mean. For example, under purely Eulerian reasoning, one might be tempted to set 
9i/u2(xi,X2,t) d 3 x 2 = so that the mean of the fluctuations would remain constant in time at a 
given point xj.. On the other hand, if the fluctuations are swept along by the mean flow, then one 
must write the Eulerian version of the Lagrangian time derivative as 

— + ui- Vij /u 2 ( Xl ,x 2 ,t)d 3 X2 = 0, (4) 

so that the mean of the fluctuation velocity remains constant along the Lagrangian mean paths of 
the ui-flow. The above relation is the second main ingredient of our approach. While the basic 
assumption of a mean flow sweeping fluctuations in a Lagrangian sense is taken as an underlying 
hypothesis, the transport equation for the mean of the fluctuation will arise as a consequence of 
our treatment. 

2.3 The kinematic model 

Applying Lemma 1 and the KSA directly to Euler's fluid equation, and setting terms in x 2 sepa- 
rately to zero yields the following equations of motion for ui and u 2 , 

(5) 

V 2 • u 2 = , (6) 

in which the pressures p\ and p 2 are determined by preservation of incompressibility of the velocities 
u 2 and u 2 , respectively. As we shall see later, the above equations also follow naturally by applying 
Lemma 1 and the KSA to Hamilton's principle. The first equation (5) yields Euler's fluid equation 
for the mean velocity Ui, while the fluctuation velocity u 2 is transported by the term (ui ■ Vi)u 2 
and also undergoes its own nonlinear Euler dynamics, given by the term (u 2 • V 2 )u 2 in the second 
equation (6). 

Equation (6) preserves the mean obtained by integration over d 3 x 2 . Moreover, integration 
of equation (6) over d 3 xi produces an averaged version of Euler's equation in the form with 

(•>i = J(-)d 8 *i 

d 

— (u 2 )i + V 2 • (u 2 ® u 2 )i = -V 2 (p 2 )i . 

Consequently, ^ J" (112)1 d 3 x 2 = 0, which is the statement of preservation of the zero-mean condition 
for the Eulerian fluctuation velocity u 2 (xi,x 2 ,t). 

Kinetic energy conservation. The total kinetic energy of a fluid flow is the L 2 norm of its 
velocity ||u|| 2 2 on the domain of flow, which decomposes according to Lemma 1 into 

Nil* = / |u 1 (x) + u 2 (x,x/e)| 2 d 3 x. (7) 

One takes the average over the rapid variations by integrating over x/e =: X2 while holding x =: xi 
fixed. This yields, 

u||£a) 2 = / |ui(xi)| 2 d 3 xi + — / |u 2 (x 1 ,x 2 )| 2 d 3 xid 3 x 2 = ||ui ||£ 2 + (||u 2 || 2 : 2) 2 , (8) 



d 

— ui + (ui • Vi)ui = -ViPi(xi) , 


Vi ■ Ui = 


d 

— u 2 + (ui ■ Vi)u 2 + (u 2 • V 2 )u 2 = 


-V 2 P2(Xl,X 2 ) 



where Vi = 1 is the volume of the domain X2 G [0, l] 3 and the zero-mean relation J U2(xi, X2) d 3 X2 = 
has been used. Thus, the mean total fluid kinetic energy decomposes into the sum of the square 
of the L 2 metric of the velocity U! and the mean-square L 2 metric of the velocity U2. 

Remark 2 The two kinetic energy norms in (7) are conserved separately by the system (5)-(6). 

Remark 3 (Relation to the alpha-models) 

The treatment so far mimics the treatment of Lagrangian averages of the WKB decomposition in 
[GH1996]. That approach led to the Lagrangian Averaged Navier-Stokes alpha model of turbulence 
[CFHOTW1998, CFHOTW1999] in which Taylor's hypothesis [Ho2005] was invoked as a closure, 
by imposing that small excitations ka > 1 evolve by being swept by the larger scales ka < 1, under 
which the nonlinearity of the smaller scales is ignored. Here, that assumption has been relaxed as 
in equation (6), to allow for the smaller scales to evolve under their own nonlinearity, relative to 
the motion of the larger scales, whose flow trajectories are treated as Lagrangian coordinates for 
the smaller scales. In particular, this means that averaging by integrating over x 2 while holding 
X! fixed may be viewed as Lagrangian averaging in this situation. 

The decoupled equations (5)-(6) comprise a simple non-interaction representation of two-scale 
Euler equations. Namely, equation (5) has reduced to Euler's equation for the velocity U! of the 
resolved scale motion (big whirls). And the velocity U2 for the subgrid scale motion (little whirls) 
evolving in equation (6) is governed by Euler's fluid equation, expressed in the moving frame of 
the Ui-flow, viewed as a scalar transformation applied to the xi-dependence of the fluctuation 
velocity 112. To describe this situation, we say that the fluctuations are swept, not slaved by the 
large scales. However, the KSA would not be enough of a nonlinear basis to describe turbulence, 
because there is not yet any back-reaction from the fluctuations to the mean motion. Having set 
up this kinematic framework, the remainder of the paper deals with modelling further dynamical 
interactions between the mean flow and the fluctuations. 



3 Formulation of the two-scale model 



3.1 Hamilton's principle for two scales of motion 

We shall derive a two-scale model of turbulence dynamics by using the relabeling symmetry in 
Hamilton's principle for ideal fluids 

5 r L(V,^,r7,77)dt = 
where we shall keep in mind the ^||u||^ 2 ) kinetic-energy form in equation (8), 

l r . ,^2 



L(il>,il>,Ti,rj) 



</>(xi 



(0). 



d 3 *i 0) 



VxD 1 * 2 



(0). 



d 3 4V*i 0) . 



Upon using the relabeling symmetry, this particular Lagrangian becomes 



L(^, -0, 77, 77) =£(ip,iJ),r}or] 1 ] 



VKxi' 



(o> 



d 3 xi 0) 



u 2 (y(xi J ),x 2 ) 



(TX 2 (TX\ 



3J0) 



where u 2 0( x i )> x 2) : = ^( x ( 0) ) ° ^Jr^X^)- One then finds that 

£(ij), ip, u 2 ) = £(ij> o if}- 1 , u 2 o tp) 

= iy'|u 1 (x 1 )| 2 d 3 x 1 + iyy|u 2 (x 1 ,x 2 )| 2 d 3 x 2 d 3 x 1 (9) 

as in (8), where ui(xi) = ^o^fxi) and u 2 (x 1 ,x 2 ) = u 2 (?/>(xi),x 2 ). 

At this point, one may take variations and apply Hamilton's principle in its general form 

5 [' £(ui,u 2 )dt = 0, 
Jtl 

for an arbitrary Lagrangian, £(ui,u 2 ) arising from relabeling symmetry arguments as above. The 
variations 8xii and 5u 2 are computed in Appendix A. 2 and they produce the following dynamics: 

d S£ „ 5£ ,„ ^ W /\ lT « 



— + Ul • Vi— + (Vi Ul ) T ■ 1 —+ (Viu 2 ) T • — - d 3 x 2 = -Vipi (10) 

Ui OUi OUi J ou 2 



«ss + «? r^J + U2 ' V2 ^ + (V2U2) ' ^ = " V2P2 (11) 

where we denote (Vu) T ■ v = VjVu^ for a vector u and a co-vector v. Equivalently, upon defining 

/ n M , , 5£ D d „ 

mill := ^— , m 2 (xi,ir 2 ) := -— , ——:= — + Ul ■ Vi , 
oui ou 2 Dti at 

we may rewrite equations (10)— (11) equivalently as 

— mi + (V lUl ) T • mi + y(Viu 2 ) T -m 2 d 3 x 2 = -ViPi (12) 

V v ' 

Div i (Reynolds stress) 
m 2 - u 2 xcurl 2 m 2 = -V 2 p 2 (13) 



Dti 



Nonlinearity 



with diviUi = and div 2 u 2 = 0, by construction. Upon specializing to the (||u||^ 2 } 2 averaged 
Lagrangian (9), one finds mi = u x and m 2 = u 2 thereby recovering equations (5) and (6). In this 
case, the Reynolds stress term in equation (12) becomes a gradient Vi / ||u 2 | 2 o? 3 x 2 , which may be 
absorbed into the pressure gradient. 

Remark 4 (Momentum conservation laws) 

Equations (12)— (13) imply the following conservation law for the advection of the total fluctu- 
ation momentum at a point xi by the mean flow, 

-^- /"m 2 (xi,x 2 )d 3 x 2 = 0. (14) 

For the kinetic energy in (9) equation (14) recovers the formula (4) for the preservation of the mean. 
In general, this momentum conservation law replaces the preservation of the mean in equation (4). 



Because rri2(xi,X2) := 5d./5\i2 is a variational derivative of an xi-translation invariant La- 
grangian, it is always possible to differentiate by parts and rewrite the last term in (fO) as the 
divergence of a stress tensor, 

(V!U 2 ) T -— d 3 x 2 = Vi-Ri. 
ou 2 

Consequently, equation (10) will also conserve the total xi-momentum, 

d 



dt J 5ui 



d d xi = . 



Likewise, equation (14) may be interpreted as a conservation law for the total fluctuation momen- 
tum 

-d 3 x 2 d 3 xi = , 



dt J J 5u 2 
whose preservation arises from Noether's theorem by X2-translation invariance. 

System (12)— (13) conserves energy and momentum, and is a Lie-Poisson Hamiltonian system, 
whose Hamiltonian function and Lie-Poisson bracket are given in Appendix A. 4. 

3.2 Evolution of fluctuation labels 

In order to compare our results with previous literature, it may be useful to compute the evolution 
equation for the fluctuation label 'Hx 1 { K 2 ) an d compare it with equation (1). This task can be 
easily accomplished by substituting the Lagrange-to-Euler map: 

^(x 1 ,X2,t) = ||^(xf\xf,0)5(x 2 -^ (x ,o, <) (xf,t))5(x 1 -^(xfU)) (15) 

into equation (11), or equivalently into equation (13). Upon denoting m 2 = 5£/5u 2 for simplicity, 

we may pair equation (13) with a divergence- less test function and obtain the following equation 

of motion 

a 

-qI»7xi( x 2 0) .*) + ui(xi,t) ■ Vir7 (x^ 0) ,t) = u 2 (r7 (x5 0) ,t), Xl ,t) (16) 

ot s « ' > v ' 

transport along the mean flow micromotion of fluctuations 

where we recall Ui = xj> o ip~ and the notation is such that xi = ^(x^ , t) and ?7 Xl (x 2 , 0) = x 2 ■ 
At this point, the comparison with equation (1) is immediate. Indeed, one readily concludes that 
the important difference between (1) and (16) is that the second allows for more freedom in the fluc- 
tuation labels, whose micromotion is encoded in the fluctuation velocity term u 2 (t; Xi (x2 , t), xi, t) 
on the right hand side. Dropping this micromotion term yields precisely (1), which assumes that 
the internal structure associated to fluctuations is completely frozen into the mean flow, moving 
with velocity u x . Actually, one can say that equation (16) embodies the essence of the present 
theory by showing explicitly how this differs from previous works in this subject. Notice that the 
equation (16) is totally equivalent to (3), although (16) is more suggestive since the right hand side 
provides a direct link with the Eulerian velocity 112 appearing in the equations of motion (10)— (11). 



4 Conservation of total SGS circulation and helicity 

4.1 Kelvin-Noether circulation theorem 

As shown in the previous section, the evolution of the fluid momentum 5£/5u 2 in terms of its 
initial value is reconstructed by using the Lagrange-to-Euler map (15). Upon dropping explicit 
time dependence for convenience, the latter can also be written by pulling back the momentum 
evolution as follows: 

xi°\xf ,0) = ||^( Xl ,x 2 ,t) *(xf ^^ 



Su 2 

so that Noether's theorem reads 

d 
di 



||^-( Xl ,x 2 ,i) sL 2 0) -V^MJ 5(xS 0) - V- 1 (x 1 ))d 3 x 2 d\ 1 = 0. 



where expanding the time derivative yields the equation of motion for 8£/8u2- 
At this point, one can fix a loop 70 and take the integral 



A / ilf x (0) x (0) 0)-dx (0) 
dt/ 70 5u 2 (Xl ' X2 ' U) dX2 

= ±<f (//^(*i>***) ^(^-^(o))^)) 5(x^-^ 1 (x 1 ))d 3 x 2 d 3 x 1 )-dxf 

to obtain the following Kelvin-Noether theorem. 

Theorem 2 (Kelvin fluctuation-circulation density) At each point x 1; consider a loop 7(112) 
moving with the fluctuation velocity u 2 (x!,x 2 ) = f] Xl or7~ i 1 (x 2 ). Then, the following transport 
dynamics holds: 

_— = where /C(x x , t) := ® - — (x x ,x 2 )-dx 2 . (17) 

Dtl / 7 ( U2 ) du 2 

The proof can be found in Appendix A. 3. 

Remark 5 Spatial integration over xi in (17) yields conservation of the total fluctuation-circulation 

d 



, /C d 3 xi = 
dt 



upon using Vi • Ui = 0. 



4.2 Vorticity and helicity density 

Applying the corresponding curl operations to equations (12) and (13) yields, upon defining 

uji = Vi x mi and o; 2 = V 2 x m 2 , 
the following vorticity dynamics 

— wi - Vi x (ui x Wl ) + ViX / (Viu 2 ) T • m 2 d 3 x 2 = (19) 

ul> 2 - V 2 x (u 2 xw 2 ) = 0. (20) 



Dt x 
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Consequently, the fluctuation helicity density 

H(xi,t) = / m 2 ■ curl m 2 d 3 x 2 (21) 

satisfies the transport equation 

This calculation implies the following. 

Theorem 3 The total fluctuation helicity is preserved, 

Proof. The proof proceeds as a direct calculation, 

j.\u d 3 xi = J -rrU d 3 xi = - J Vi(Tiui) d 3 xi = 
for periodic boundary conditions. ■ 

5 Resolved- scale models based on SGS advection 

5.1 A resolved-scale model that is structurally similar to MHD 

The advection laws (17) and (22) for the total SGS circulation K and helicity "H, respectively, 
suggest replacing the SGS degrees of freedom with these Lagrangian-averaged quantities obtained 
from X2-integration. Hamilton's principle then takes the form 

Jti 
This variational principle results in the following equations of motion for rri! := 8£/8ux, H and /C, 

^i + (V lUl ) T • Hi, + ^V,W + ^V,K = -V lPl (23) 

^ = 0. ^ = 0- (24) 

Dti Dt x 

An illustrative special case can be found by choosing the Lagrangian 



£(u u H,K) = \ /luiCxOI 2 - h^xU x V!/C| 2 d 3 



x x . (25) 



Upon defining B := Vi"Hx Vi/C, the corresponding equations (23) and (24) arising from Hamilton's 
principle for this Lagrangian are 

^i + Bxcurl 1 B = -V lPl , Vi-ui = 0, (26) 

L)tx 

— = curl 1 (u 1 xB), Vi-B = 0. (27) 

11 



These are identical to the equations of ideal incompressible magnetohydrodynamics (MHD). How- 
ever, in these equations, the quantity B is the cross product of the xi-gradients of the X2-integrated 
SGS helicity % and circulation densities /C, defined in equations (21) and (17), respectively. The 
equations in (27) imply the relation 

D 



B • ViUi , so 



B- Vi 



0, (28) 



Dti IDti 

where [ • , • ] is the commutator of divergence-free vector fields. Consequently, one finds the relation 

^5 = _*L(b • V lUl ) = - B ■ Vi (V x ( Pl + I|B| 2 ) - B ■ V x b) . (29) 

Remark 6 The two formulas (28) and (29) correspond to Ertel's theorem [Erl942] and Ohkitani's 
relation [Ohl993] for Euler's equations with vorticity w = V x u, respectively, 



D 



and -^ = -(w-Vu) = -wVVp. (30) 



Remark 7 Equation (29) may be rearranged into the form of a nonlinear wave equation 

^2 - ( B ■ Vx) 2 B = - B ■ ViVi (Pi + ||B| 2 ) . (31) 

We have shown that the SGS advection laws (17) and (22) for the integrated SGS circulation 
K. and helicity % impart a certain elasticity to the resolved scale equations, which is completely 
analogous to Alfven waves in MHD. The magnitude |B| is not preserved by this flow. In fact, 
equation (28) implies 

« = B .(V 1 u 1 ).B (32) 

for the evolution of the magnitude. 

Remark 8 If the Lagrangian in (25) were actually relevant in turbulence, the result (31) would 
have vast implications for the spatiotemporal properties of turbulence. In particular, the cross 
product of resolved gradients Vi/C and V{H of the integrated SGS circulation and helicity would 
propagate as a wave. This is a degree of freedom not usually considered in turbulence models. 
However, we have treated the Lagrangian (25) here just for illustration. We do not expect it to 
actually describe turbulence. We have discussed it only to show that an energetic dependence on 
the resolved gradients of the integrated advected properties of the subgrid scales could lead to 
interesting dynamical behaviour in the resolved scales, by imparting a type of potential energy 
that would lead to a type of non- Newtonian behaviour. The concept of non-Newtonian properties 
of turbulence has a long history of investigation, going back at least to [R11957]. More recently, the 
related concept of nonlinear dispersion in turbulence has also been studied intensely, see [FHT2001, 
FHT2002]. However, as far as we are aware, the closest analog in the previous literature of this 
concept of elasticity of integrated SGS properties is due to [MP1997], who derived similar equations 
from a different approach, based on stretching of SGS vortices. A derivation of a class of stretched- 
vortex SGS models similar to those discussed in [MP1997] will be provided using a variational 
principle in the next section. 

Remark 9 The two formulas (28) and (29) correspond to Ertel's theorem and Ohkitani's relation 
[Ohl993] for Euler's equations with vorticity u> = V x u, respectively, 



D 



and — ^ = — («-Vu) = -w VVp. (33) 

L)t z Dt 
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5.2 Convection in the stretched SGS vortex model 2 of [MP1997] 

In this section, we introduce the constrained kinetic energy Lagrangian 

£\ Ull 2_ p{D _ l) _ qilBl 2_ 1)d 3 x ^ (34) 

with B := Vi"H x Vi/C and Dd 3 Xi the preserved volume element. The Lagrange multiplier p (the 
pressure) imposes volume preservation. The Lagrange multiplier q imposes |B| 2 = 1, so that, as in 
[MP 1997], we may think of the subgrid scale order parameter as a vortex filament whose strength 
is constant, but whose direction varies. The Euler-Poincare motion equation for this Lagrangian 
is obtained as in [HMR1998] and found to be 

Dux „ „ 2 „ n ( r , B<g>B 



Vip - gVilBl 2 - Vi • 2q (Id - -^r (35) 



MP97 SGS stress form 

The motion equation (35) conserves the kinetic energy | J |ui| 2 d 3 x\. This means that, unlike 
most SGS models of turbulence, the subgrid scale vortices in the present model do not dissipate 
resolved-scale kinetic energy. Instead, all dissipation of kinetic energy at the resolved scales must 
be modelled separately; for example, as a viscous term that could be added later at the xi-scale. 
Models of dissipation are not part of our discussion here. 

The last expression in equation (35) has the same form as the stress term in equation (25) of 
[MP 1997], although the corresponding expression there does not transform properly as a tensor. 
Instead, the expression for the stress tensor there transforms under change of variables as an array 
of scalars. This means it neglects the convection of the subgrid vortices by the resolved field, as 
discussed explicitly in [MP1997]. It also means that the stretched SGS vortex model 2 of [MP1997] 
is not variational. The present variational model also contains auxiliary equations for preservation 
of volume and advection of B by the resolved field, with constant magnitude, 

— = -divi(Dui), — = curli (ui x B) , |B| 2 = 1. (36) 

at at 

The Lagrange multipliers p and q obey a system of linear equations found by preservation of initial 
conditions D = 1 and |B| 2 = 1, which require, respectively 

B • ( Villi) ■ B = with Vi-B = 0. (37) 



1 D D- 


= Vl • Ul = 


, 1^|B| 2 


DDt x 


2 Dt x 


Summary 







(i) When q is interpreted spectrally as 2q(t, Xi) =: /w tx ,E{k)dk for a cut-off wavenumber 
k c (t, Xi), equations (35)-(37) represent a variant of the MP97 model 2 in which we have 
restored the convection of the subgrid vortices by the resolved field that was neglected in 
[MP1997]. 

(ii) When 2q = 1, this model reduces to ideal incompressible MHD, in which case D\Q\ 2 / Dti ^ 0, 
by equation (32). 
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Ertel Theorem and Ohkitani relation for the present variant of MP97 model 2. The 

present variant of the MP97 model 2 motion equation (35) may be written as 

^Hi = -Vi(p + 2g) + V 1 -2gB®B=:F ! (38) 

where preservation of Vi • ui = and B • ViUi ■ B = determines p and q, and 

<9B 

— = curli (ui x B) , |B| 2 = 1 with Vi • B = . (39) 

The Ertel Theorem and Ohkitani relations for this variant of MP97 model 2 are then 



D 



, B Vx 



, for 2— = B ■ ViUi and ^ = -£-(B ■ V lUl ) = B ViF (40) 
L)t\ LJX\ L)t\ 



.Dti 

Together, the Ertel and Ohkitani relations conveniently deliver 



D 



(B • ViUi • B) = B • ViF B + IB ViUi 



2 



Dtx 

which must vanish in order to preserve the constraint that B 2 = 1. 

The equation system for Lagrange multipliers p and q. Preservation of Vi ■ ui = and 

B ■ ViUi ■ B = determines the Lagrange multipliers p and q from the following system 

d .- , ,_ |2 



— (V 1 -u 1 ) = -|V 1 u 1 | 2 + V 1 -F = 



D 



(B ■ ViUi B) = B ViF B + |B ViUi 



2 



(41) 



Dt x 

in which F := -Vi(p + 2q) + Vi ■ (2gB <g) B) . 

That is, the force F for this variational version of MP97 model 2 depends linearly on p and q as 
in (38). An additional g-term enters the boundary conditions for this system. Namely, the normal 
component of the force F must vanish on a fixed, flat boundary; that is, F • h = 0. 

Remark 10 See the papers [MP 1997, VPC2000] for discussions of numerical implementations 
of the stretched-vortex subgrid-stress model, as well as discussions of its applications to forced 
and decaying turbulence, and studies of its numerical complexity relative to standard Large Eddy 
Simulation (LES) models such as the Smagarinsky model. The present variational version of the 
stretched-vortex subgrid-stress model restores convection of the subgrid vortices by the resolved 
field. However, it has the additional complexity that the system of equations (41) for p and q must 
also be solved at each time step. The solvability of this system has not been studied yet and such 
a study would be beyond the scope of the present work. 

6 SGS fluid model with feedback between the scales 

6.1 Energy coupling of fluctuations and stress tensor 

In this section, we shall return to the multiscale description and consider the particular Lagrangian 

£(u 1 ,u 2 ) = iy'|u 1 (x 1 )| 2 d 3 x 1 + iyy(|u 2 (x 1 ,x 2 )| 2 + ^ 2 |V lU2 | 2 )d 3 x 2 d 3 Xl , (42) 

Coupling term 
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where ct\ is an appropriately chosen coupling constant with the dimensions of length. In this case, 
the mean-fluctuation interaction energy is the L 2 norm of the xi-gradient of the slowly varying 
envelope of the fluctuation velocity u 2 (x!,x 2 ). The momenta mi and m 2 in equations (12)— (13) 
are given in terms of the velocities ux and u 2 for this Lagrangian as 

rri! = u x and m 2 = (1 — a 2 Ai)u 2 , hence u 2 = (1 — a^A 1 ) _1 m 2 , 

where Ai is the Laplacian operator in the xi coordinates, so that (1 — a^Ai)" 1 is a smoothing 
operator. Consequently, the coupling term in (42) introduces a Reynolds stress term that will 
remain and affect the evolution, as follows. 



D 
Dt 



- Ul - o^Vi • / (ViuJ • Viu 2 - Id|Viu 2 | 2 ) d 3 x 2 = -Vipi (43) 

1 ' 1 . ' 

(Divergence of Reynolds stress, diviR iy ) 

'u 2 - c^A!U 2 ) - u 2 x curl 2 (u 2 - c^Aii^) = -V 2 p 2 . (44) 



Dti 



(Nonlinear convection) 
Upon introducing the subgrid Reynolds stress notation R 2 , the latter equation becomes 

D 



Dt x 



[u 2 - a?A lU2 ) + V 2 • (r 2 - Id p 2 ) = (45) 



(Subgrid Reynolds stress) 



where the term 



u 2 x curl 2 (u 2 - c^Aiu 2 ) = V 2 • R 2 - u 2 (V 2 • u 2 ) (46) 



has zero mean for V 2 • u 2 = and the transformation to V 2 • R 2 stress-divergence form in (45) 
represents momentum conservation, which arises from Noether's theorem for the x 2 -translation 
invariant Lagrangian in (42). 

1. Reynolds stress. The only channel of feedback from small scales to the larger ones arises 
from the Reynolds stress term in equation (43). The Reynolds stress tensor 

ViuJ • Viu 2 d 3 x 2 with components RV —a 2 — — \ — — § d 3 x 2 (47) 

J ox* ox\ 

in equation (43) is remarkably similar to the tensor diffusivity stress tensor for large-eddy 
simulation (LES) introduced in [Lel974] and discussed in detail in [WWVJ2001]. It ex- 
changes energy between the smaller scales and the larger scales along the L 2 mean primary 
stretching direction of Viu 2 , and the direction of this exchange depends on the gradient of 
the slowly varying envelope of the fluctuation velocity. There is, however, one very impor- 
tant difference between the Reynolds stress in (47) and the corresponding term in an LES 
model: our model does not yet take viscous effects into account. Indeed the system (43)- 
(44) conserves the kinetic energy expressed by the Lagrangian in (42), and even comprises a 
Lie-Poisson Hamiltonian system, as discussed in Appendix A. 4. 

The Reynolds stress term in our system represents backscatter due to nonlinear dispersion, 
not diffusion. Thus, as with the Euler-alpha model [HMR1998, FHT2002] the higher order 
terms that make the solutions more regular do so without introducing dissipation. 
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2. SGS circulation and helicity advection. Note that the present model retains the same 
advection laws (17) and (22) for the total SGS circulation and helicity as in Section 4. 
Therefore, the corresponding resolved-scales exist and may be studied for this choice of the 
multiscale Lagrangian, as well. 

3. Taylor hypothesis. The Taylor hypothesis for the small scales assumes that their nonlin- 
earity does not affect their evolution; so that they are passively swept by the larger scales. 
Applying the Taylor hypothesis would neglect the term marked Nonlinear convection in equa- 
tion (44), or Subgrid Reynolds stress in equation (45). Even if it eliminated the nonlinear 
convection effects of the fluctuation velocity, the Taylor hypothesis would not eliminate the 
effects of the back-reaction of the small scales due to their sweeping by the large scales, which 
still remain in the term labelled as Reynolds stress in equation (43). 

4. Regularity of the equations. Equation (44) is reminiscent of the Navier-Stokes-alpha 
model [FHT2002], but with smoothing applied to the envelope of the fluctuations in a frame 
moving with the mean flow. In fact, equation (44) is essentially a two-scale version of the 
Euler-alpha model [HMR1998] in a frame moving with the xi-fiow. 

The Lagrangian (42) could be modified further to allow introduction of norms that would 
be strong enough to ensure long-time existence of its corresponding solutions, even in the 
absence of viscosity. In the presence of viscosity, the resulting equations are: 

-Ui + diviRi = -Vipi + V1A1U1 , (48) 



— — (u 2 - a^Aiu 2 ) + div 2 R 2 = -V 2 p 2 + z/ 2 A 2 (u 2 - a^Aiu 2 ) , (49) 

L)t\ 

with R 2 and Ri given in (46) and (47), respectively. The flows at both scales are incompress- 
ible, so the pressures pi and p 2 in equations (48) and (49) are determined from preservation 
of diviUi = and div 2 u 2 = 0, respectively. The viscosity has been introduced ad hoc here, as 
the diffusion of momentum at the each scale. Equations (48)-(49) may admit global strong 
solutions, as occurs for the Navier-Stokes-alpha model investigated in [FHT2002]. However, 
their analysis is beyond the scope of the present paper and will be left to the future. 

5. Numerical complexity and implementation of multiscale models. It is possible that 
one may be able to reduce the computational complexity of these equations and thereby 
accelerate their computation by introducing a type of optimal sampling that would select 
those smaller scales that make the largest contributions to the Reynolds stress term. The 
multiscale analysis in the present case couples the resolved scale solution with a subgrid cell 
problem for vortex filament evolution at each point of the resolved scale grid. The computa- 
tional cost for this coupled system of equations could be quite expensive, although there are 
some alternatives that might be used to lessen the cost. For example, an adaptive scheme 
recently has been developed to reduce complexity and speed up the computation in a related 
case. See [HYR2008] for a discussion of a promising approach for reducing the complexity of 
such multiscale computations. See also [EH2009] for an approach to multiscale computations 
using Finite Element methods. Yet another possibility might be to adapt the heterogeneous 
multiscale method [HMM] that has already been developed for complex fluids to the present 
case. The difference is that in the present case the subgrid scale dynamics governs an infinite- 
dimensional vortex filament, rather than a finite-dimensional order parameter, as occurs in 
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liquid crystals and magnetic fluids. In any case, the numerical simulations of the equations for 
the present multiscale model and the development of numerical algorithms for their solution 
must be left to the future. 

7 Conclusions 

Using a simple decomposition argument that applies to all periodic functions, the present approach 
has sought to transfer ideas associated with convection of microstructure in complex fluids into 
the context of turbulence modeling. The microstructure was interpreted as an order parameter in 
the same fashion as spin is regarded for ferromagnetic fluids. This interpretation was suggested 
by previous work in which turbulence microstructure had been assumed to be transported by 
the mean flow. However, in this paper, the intrinsic nonlinear features of microstructure were 
also considered so that microstructure underwent its own nonlinear evolution in the frame of the 
mean flow. In short, the kinematic sweeping ansatz (KSA) assumed that the mean flow serves 
as a Lagrangian frame of motion for the fluctuation dynamics. This seems to be an effective 
approximation, in general, and it was contrasted in section 3.2 with the exact formula. Other 
assumptions about convection of microstructure would have been possible. However, the KSA 
has the advantage of possessing a purely Eulerian description, and thereby avoiding the necessity 
of computing the Lagrangian back-to-labels map. In addition, the geometric framework that 
was previously developed for complex fluids could be naturally transferred to turbulence, thereby 
leading, in the present case, to pure transport dynamics for the fluctuation-circulation and the 
fluctuation-helicity densities. 

The resulting model arising from KSA has also the advantage that additional, even finer, scales 
could also be incorporated by simple extension of the geometric features. In the resulting iterated 
hierarchy of smaller and smaller scales, the flow associated to a certain scale serves as a Lagrangian 
coordinate for the flow of the next finer scale, and so on. 

Even in the simplest 2-scale version treated here, the KSA yields dynamical equations that we 
hope will be suitable for modelling purposes. For example, when the lengthscale a is inserted in 
Hamilton's principle for smoothness requirements, the resulting model produces Reynold's stress 
tensor in a way that is quite reminiscent of the diffusivity stress tensor in LES simulations. 

Various features of the present model still remain to be addressed. For example, vortex dy- 
namics at the smaller scales could be an interesting future direction for research, in both 2D and 
3D. Indeed, although the ux equation may lose the vortex filament solution, the u 2 has a vor- 
tex filament solution at each Eulerian point x 1; see Appendix A. 5. Hence, projecting onto the 
plane would yield a proliferation of point vortices at each point in physical space. In addition, 
the development of numerical integrators that respect the geometric framework here would be a 
very important advance, both in new mathematics and toward evaluating the present model in 
applications. For discussions of recent progress in this direction see [GBPav2012]. Finally, we hope 
that the methods of [HMM] might also be profitably applied to the present multiscale model of 
ideal fluid motion. 
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A Appendix 

A.l Proof of Lemma 1 

Consider the Fourier series expansion of a periodic function u(x) on the unit cube [0, l] 3 , 

u(x) = J2 u(k)e ik - x , k e Z 3 
keza 



?:', 



in which u(k) are the Fourier coefficients. Let < 1/N < 1 to be a reference wavelength, where 
N is an integer, and one denotes 

A N = {k e Z 3 ; |k| < N} and A'^ = Z 3 \Atv 

Decompose the function u into two additive parts as follows: 

u(x) = ui(x) + u 2 (x), 

where 

Ul (x) = ^ u(k)e ik ' x , and u 2 (x) = ^ u(k)e*- x . 

keAiv kGA^ 

Rewrite k = kx + Ak 2 where kx and k 2 take integer values in Z 3 with k x e A^v, and compute 

u 2 = J2 u(k)e lk ' x 

kGA^ 

J2 u(k! + Ak 2 )e i(kl+iVk2) - x 

(ki+Afk 2 )GA' JV 



J2 ( Yl u(ki + Ak 2 )e iklX ) 



ik 2 -iVx 



=: ^u 2 (k 2 ,x)e Mx 

=: u 2 (x,x/e) 
where e = 1/N and the quantity 

u 2 (k 2 ,x):= Y Q(ki + Ak 2 )e ikl - X , 

kiGAjv 

involves only Fourier modes whose wave number is less than N in magnitude. Hence, any periodic 
function uel 3 may be rewritten as 

u(x) = ui(x) + u 2 (x,x/e) = Ui(xi) + u 2 (xi,x 2 ) 

where u 2 (xi,x 2 ) is a periodic function in x 2 with mean zero, 

u 2 (xi,x 2 )rf 3 x 2 = 0, 

since k 2 ^ 0. Note that the functions ui and u 2 depend on the choice of e = 1/N. ■ 



Remark 11 The decomposition 

u(x) = ui(x) + u 2 (x,x/e) 

= J2 u(k)e ik ' x + ^u 2 (k 2 ,x)e' k2X / e 

keAiv k 2 ^0 

may be regarded as a Fourier- series generalization of the WKB form 

u(x,x/e) = u(x) + - (a(x)e i9 W/ e + a*(x)e" ie(x)/£ ) 

This WKB form was used in [G HI 996] to develop a wave, mean flow interaction theory by apply- 
ing Lagrangian averaging in Hamilton's principle for rotating, stratified incompressible flow. The 
present work is similar in spirit to that previous work. 

A. 2 Expressions for the velocity variations 

The Euler-Poincare equations are obtained from Hamilton's principle with 

Sui =£i- ad ui £i 

<5u 2 = £ 2 — ad n2 £ 2 + <£ Ul £2 — £^2 

= 6 - ad U2 £ 2 + (ui • Vi) £ 2 - (£i ■ Vi) u 2 

where diviUi = 0, div 2 u 2 = 0, £ w v = w • Viv denotes the Lie derivative, and 

adt«6 = u 2 ■ V2C2 - 6 • V 2 u 2 , 

In particular, while the first variation follows directly by taking Sui = 5( , j/>?/>~ 1 ) and by defining 
£1 = <5?/>i/' _ , the second variation follows by the calculation below, which uses the pullback 
notation: 

= £^u 2 - u 2 ■ V 2 £ 2 + ( 2 - £ Ul i 2 + i 2 ■ V 2 u 2 
where we have denoted £ 2 = ^*[(5(ip*rj)) (i/?*ry)~ ). 

A. 3 Proof of the Kelvin-Noether theorem 

Upon denoting 

S£\ 5£ , . x ( 81 \ S£ , ( o) (o) 



-1 



x l5 x 2 ,t), — =_(xr,x^,0) 



K 8u 2 J t Su 2 \5u 2 J <5u 2 

and by using the pullback notation, we have 

£)r*<™-(&-(&-™-*(£,. 
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Thus, we can take the circulation around a fixed loop 70 in R| 



°= 5 £<*•*»•**(£;, 



^70 

dt\ J 1 \Su 2/t 



SI 
<fri 2 



*-|£ 



dt J <fri 2 



dx 2 + (ui • V] 



A <fri 2 



dx- 



where 

moves with velocity u 2 and 



7 = (^*^)°7o = ^(xi)(7o) 



7 = V>*7 = »fy-io^(xi) = 7 ?x 1 (7o) 
moves with velocity u 2 = if>*u. 2 . Thus, upon applying t/» # , we have the Kelvin circulation theorem 

— * dx 2 + (ui • Vi) * dx 2 = . ■ 

at y 7(u2 ) ou 2 J 7(U2) ^u 2 



A. 4 Hamiltonian structure 

Upon defining nij = S£/SUi, one writes the functional Legendre transformation 

/i(m b m 2 ) = / mi ■ Ui d 3 Xi + / / m 2 • u 2 d 3 x 2 d 3 Xi - £(u x , u 2 ) 



so that the equations of motion (10) and (11) read as 



<9mi Sh 



dt 5m.i 



■ VlHli 



„ Sh 

Vi- — ] -mi 

omi 



/«)*■ 



m 2 d x 2 = — V1P1 



9m 2 d f Sh \ Sh 
at dx{ \Su{ / ^m 2 



'<frn 2 / 



m 2 = -V 2 p 2 . 



Therefore, by standard methods [HMR1998], the Poisson bracket for the above system is 



if,h} 



mi ■ 



Sf Sh 



Sm.i <krii 



d 3 xi 



m 2 ■ 



Sf Sh 



+ / / m 2 ■ 



(frn 2 (frn 2 



d 3 x 2 d 3 xi 



5mi <5m 2 5m 1 (frn 2 

where we have denoted the Lie algebra brackets as 

[v, w] . = w ■ ViV - v ■ VjW 
fori = 1,2. 



(50) 
(51) 



•Vi^-^.Vi^)d 3 * 2 d 3 *i 
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A. 5 Vortex structures in the subgrid scales 

This appendix shows how equations (43)- (44) allow for vortex structures in the fluctuation dy- 
namics. In order to see how this happens, we can write (12)-(13) in terms of the vorticities 

u*i = Vi x mi , u>2 = V2 x m 2 

which yields (19)- (20) with 

/"(Viu 2 ) T • m 2 d 3 x 2 = - /(Viu; 2 ) T • <f> 2 d 3 x 2 (52) 

where <p 2 is the vector potential associated to u 2 , so that u 2 = V 2 x </> 2 - 
Notice that equation (20) possesses a vortex solution of the type 

w 2 (xi, x 2 , t) = 9 a R(xi, s, t) 5(x 2 - R(xi, s, t)) ds , 

so that a vortex filament in the fluctuation vorticity is attached to each point of the mean fluid. 

In turn, the above vortex solution can be used to reduce the level of difficulty of the equations 
(12)-(13). Indeed, replacing the above solution into (52) eliminates the integral over x 2 thereby 
yielding 



J (V^ 2 ) T • 2 d 3 x 2 = - j (0 2 ■ ViR' - R V 2 2 : ViR) ds 



X2=R(xi,s,i) 



where R' = <9 S R and A : B = AijBij. The above term shows how the subgrid vortex determines 
the strain tensor of the mean flow. On the other hand, the dynamics of the vortex filament can 
be derived by pairing equation (20) with a test vector field, thereby yielding 

(dt + Ui • Vi)R= V 2 x 02lxa=R(xi,M) = u 2(xi,R(xi,s, £),£), 
which shows how the mean velocity affects the subgrid vortex via the material time derivative. 
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